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Abstract 

We discuss the characteristic interference features of soft radiation in the 
threshold production of heavy unstable particles: soft gluon radiation in e + e~ — ► 
ti and soft photon radiation in e + e~ — > W + W~ . We show that the heavy 
particle decay width controls the interference between the emission off the final 
state particles. As a result, the radiation pattern may provide a way of measuring 
the decay width of the heavy particles. 



1 Introduction 



Heavy unstable charged particles can emit radiation both before and after they decay. 
The analysis of such radiation is a complex issue, depending sensitively on the timescale 
of the emission compared to the lifetime of the unstable particle [1]. In particular, 
the radiation pattern can be very different according to whether the radiation occurs 
predominantly during the production stage or after the particle has decayed [2]. 

There are several important examples of such effects which are directly relevant to 
present and future high-energy colliders. As a specific example, consider the produc- 
tion and decay of a ti pair in high-energy e + e~ annihilation. With a mass of at least 
91 GeV [3], the top quark can decay to a real W boson and a b quark. The width 
r t for this decay is quite large — so large that the top weak lifetime can be as short 
as strong interaction timescales. The resulting interplay between the strong and weak 
interactions of the top quark gives rise to interesting physical effects. For example, if 
top is heavier than ~ 100 GeV, then T t can be greater than the typical hadronic scale 
ji ~ 1 fm _1 and it may decay before it has time to hadronize [4-6]. In particular, ti 
resonances may never be formed. Here we are interested in the perturbative aspects of 
the strong-weak interplay: decay versus gluon bremsstrahlung. Reference [2] discussed 
soft gluon radiation in e + e~ — > ti and showed that gluons radiated in top production 
and decay can interfere, and how much they do depends on the top width. This means 
that top production and decay should not be treated separately - the gluon distribu- 
tion in top events is not what one might naively guess. Furthermore, this width effect 
might be useful; the sensitivity of the soft gluon distribution to T t suggests a way 
to measure it [2,7]. The width dependence of the gluon distribution at high collision 
energies was studied in Ref. [2]; however, it was found that the configurations with 
the most sensitivity to T t were also the least likely to occur. 

In this paper we consider on soft gluon radiation near the ti threshold. The top 
quarks are produced nearly at rest and essentially do not radiate. The width depen- 
dence is a result of interference between gluons radiated in the two decays, which 
does not play an important role at higher energies. Near the production threshold, 
the amount of interference between gluons from the b and b is controlled by the top 
width, and what matters is the size of T t relative to the gluon energy. Thus we will 
see that when the top width and the gluon energy are more or less the same order of 
magnitude, the radiation pattern is sensitive to T t . 

A second process which exhibits similar features is the emission of photon radiation 
in the process e + e~ — > W + W~ — > //'//'■ The radiation pattern of a soft photon of 
energy oo is sensitive to the W decay width for uj ~ Tw Here there is the additional 
complication of radiation off the initial state as well, but, as we shall see, this can 
easily be taken into account. 

In principle, therefore, the study of the soft gluonic and photonic radiation in ti 
and W + W~ production provides a basis for determining the decay width of the heavy 
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particle. In practice, however, there are many difficulties. The measurement of the 
radiation pattern in top quark production requires the separation and identification of 
a soft gluon jet (typically with energy u ~ 5 GeV). In the case of photons radiated in 
W + W~ production, while identification of relatively soft photons might not pose too 
many problems, the event rates are low for the anticipated luminosities of future e + e~ 
colliders. Nevertheless, we believe these issues are worth exploring for several reasons. 
First, on a theoretical level there are several features of the radiation patterns that 
show interesting interference effects which are at first sight counter-intuitive. Second, 
the 'traditional' methods of measuring the masses and decay widths by threshold 
scanning are not without their own problems. Especially for tt production at high- 
energy e + e~ linear colliders, the structure of the threshold is smeared by beam-induced 
effects, intrinsic energy spread etc. On the theoretical level, the measurement of T t 
from the shape of the cross section as a function of beam energy near threshold is a 
delicate issue. For e + e~ — > W + W~ , the threshold scan strategy requires a detailed 
calculation of higher order electroweak effects including width effects and initial state 
radiation. To our knowledge, a comprehensive one-loop calculation including the width 
effects has not yet been completed. An attempt to incorporate Yy/ into the tree- 
level formulae has been made in Ref. [8] (see also Ref. [9]). These results however 
are seriously affected by initial state radiation and other effects like the final-state 
Coulomb attraction near the W + W~ threshold [10, 11]. As a general comment, one 
might argue that the W width is already known to fairly high precision (IV = 2.15 ± 
0.11 GeV [12]) from indirect measurements using W production cross sections in pp 
colliders. Nevertheless, we believe it is important to obtain a direct measurement of 
this important Standard Model parameter as well. 

In this paper, therefore, we will focus mainly on the theoretical features of soft 
radiation in tt and W + W~ production. Some illustrative numerical results for tt 
production are presented; a more complete numerical treatment of the W + W~ case, 
particularly with reference to LEP200, is deferred to another paper [13]. We believe 
that our conclusions will show that a more detailed experimental study (event rates, 
detector capabilities, etc.) is certainly warranted. The remainder of the paper is orga- 
nized as follows. In Sections 2 and 3 we discuss the radiation pattern, for gluons in tt 
production and photons in W + W~ production respectively, near threshold in detail. 
In Section 4 we present numerical results for top production and discuss prospects for 
measuring r t . We conclude in Section 5. Appendices contain a semi-classical deriva- 
tion of the radiation pattern and further details of the calculation of the distributions. 

2 Soft radiation pattern 

We are interested in emission of a gluon in the process e + e~ — > tt — > WWbb and of a 
photon in the process e + e~ — > W + W~ — > //'//'■ Although the analysis of the final 
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state radiation in both processes is very similar, the latter process is complicated by 
the additional contributions from initial state radiation. To begin with, therefore, we 
discuss the tt case, and extend the analysis to W + W~ production in the next Section. 

The general result for soft gluon radiation in e + e~ — > tt — > WWbb was presented in 
reference [2] (see also [7]). Here we focus on the particular case of radiation close to the 
tt threshold. There are two advantages in this. First, the production cross section is 
largest just above threshold. Second, near threshold the top quarks are almost at rest 
and only the b and b radiate. While it is not obvious that the top quark width should 
enter at all if only the b-quarks radiate, we can understand the its role as follows. 
Consider two cases of gluon radiation from a bb pair. If the quarks could radiate 
independently, with no interference, the gluon distribution would be proportional to 

-o -o -o u i sin2 ^ v % sin2 ^2 

(1 — fiCOSfli) 2 (1 — 1> 2 COS0 2 ) 

where Vi is the velocity of the b (b), and 0, is the angle between the b (b) and the gluon 
for i — 1 (2). (Note that in what follows we will make a distinction between v\ and 
v 2 although in practice, for the case of e + e _ — > tt — > W + W~bb in the centre-of-mass 
frame, we always have v± — v 2 .) In the other extreme, with interference we have 
coherent emission, and the gluon distribution looks like 

J = 



viv 2 {cos 81 COS 6 2 - COS #12) . . 

(1 — V 1 COS0i)(l — f 2 COS0 2 ) 



where Q\ 2 is the angle between the b and the b. This is just the familiar antenna pattern 
for emission from a quark- ant iquark pair; the interference is 2J. Note that these 
patterns can be quite different and the interference can be constructive or destructive. 
In particular, and as we shall discuss in more detail later (and see Appendix E), 
coherent emission exhibits angular ordering behavior, i.e. if we integrate over the 
azimuthal angle about the direction of quark 1, all radiation from quark 1 is suppressed 
for angles Q\ such that [14,15] 

v 2 cos 0i2 < cos 9i < cos 9 C ~ vi . (2.3) 

The factor vi on the right-hand-side takes account of the screening effects due to the 
mass of quark 1 (dead-cone). 

Now recall that in the case of interest the b and b are produced by the decays 
of the t and i. If the top lifetime is very short compared to the characteristic time 
for emitting a gluon of energy u (i.e. T t ^> u), the b and b are produced nearly 
instantaneously, and we expect coherent gluon emission, with the gluon distribution 
determined by lZ cohcI .. If however the top lifetime is very long (i.e. T t <C cu), the b and 
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b are produced at very different times and thus will radiate independently; the gluon 
distribution is then given by 7^ iIldcp .. The ratio T t /uj of the top width to the gluon 
energy controls the amount of interference. The full distribution we show below was 
derived from standard Feynman diagram techniques in Ref . [2] . But in fact it can be 
derived from simple semi-classical wave arguments - see Appendix A. 

Following reference [2], the radiation pattern can be presented as a probability 
density, normalized to the lowest order cross section: 

, / , duo dQ CpOt s . . , n , s 

dN = 1 a da q = — — j\f . 2.4 

uj An it 

where Cp = 4/3 is the QCD colour factor. Near threshold, we have 

J\f EE (1 - • "^indcp. + • ^cohcr. 

= TZ indcp . + 2 X (uo) ■ J , (2.5) 



where lZ indcp . and TZ cohov . are given in eqns.|2~T| and ^]2] (with v\ = v 2 = u&) and we have 
introduced the profile function 

xH = ™z-5 » ( 2 - 6 ) 

The factor x, which depends on the decay width and gluon energy, determines the 
amount of interference, and as stated above we have in the limits of large and small 
width, 

Af = TL indcp , for c;>r, (2.7a) 
N = TL co ^. for w<r. (2.7b) 

Evidently, we have maximal sensitivity to T t for uj ~ T t , and this provides a possible 
basis for measuring the width. For gluon emission we must in addition impose uj > 
/i ~ 1 fm~ , in order to remain in the perturbative regime. 

Before performing a detailed numerical study of the above to investigate the ac- 
tual sensitivity to T t in different angular configurations, we make some additional 
comments. 



(i) The numerator of the J term in Q2.2|) can be written as 

cos 6i cos 9 2 — cos 9 12 = — sin 9\ sin 9 2 cos (f>i 2 (2.8) 

with 0i2 the relative azimuth between pi and p 2 with respect to the gluon direc- 
tion, k. Thus an immediate consequence of the interference is that the azimuthal 
symmetry of the radiation about the 6-quark directions is destroyed. Note also 
that the expression (|2.8|) vanishes when the direction of the gluon momentum is 
chosen to be close to that of one of the quarks, 9\ 9 12 or 9 2 9 12 . Therefore 
the ^-dependence of the radiation pattern (|2.5| ) can reveal itself only if the gluon 
emission angles are not small compared to the opening angle of the b and b . 
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It is fairly straightforward to integrate over the angle of the emitted gluon and 
obtain the total probability for the radiation of a gluon of given energy. (In 
a sense this is a formal procedure, since in practice the gluon jet will only be 
identified when separated from other final state jets and from the beam.) Full 
details are given in Appendix B - here we present only the result: 



dN 
duj 



C F a s 



71UJ 

+X 



i-xM) 

1 + r 



1 , 1 

— log T 



Vl , 1 , 1 

— + — log 

V 2 1 — V 2 



-log 

r 1 — r 



where 



1> 



Pi 



4 

■P2 



MiM 2 



'1—Viv 2 cos 9 12 ) 
(l-vf)(l-vl) 



(2.9) 



(2.10) 



Note the logarithmic collinear singularities which dominate the integrated dis- 
tribution in the ultra-relativistic limit t>, — > 1. In fact for (1 — Vj) <C 1 and for 
612 values not particularly close to 0, we have (Appendix B) 



dN 



C F a s 



-RUJ 



In 



1- Vl 



+ In 



I—V2 



4 + 2 X (oj) 



In 



cos 6 



12 



1 



(2.11) 



Note that the second tu-dependent term in fl2.11|) enhances or depletes the radia- 
tion according to whether 612 is larger or smaller than @ crit w 75° (Appendix B). 



3 Soft photon radiation in e + e — > W + W 

Before presenting our numerical results for the ti case we discuss in this section the 
extension of the above analysis to soft photon radiation near threshold in the process 
e + e _ — > W + W~ — > ff'ff'- Once again the radiation pattern includes contributions 
from the WW production and decay antennae, together with interferences between 
them. As was demonstrated in Ref. [16], the radiation at the production stage (the 
WW antenna) is given by the classical current expression, as for the ti antenna [2], 
irrespective of the choice of gauge. Apart from overall couplings, colour factors and 
charges, there are only two main differences. First, there are additional contributions 
from initial state radiation. These pose no particular problems in practice as long 
as the final state particles, including the photon, are kept well away from the beam 
direction. As shown in reference [17], when all the final state particles are exactly 
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transverse to the beam direction, the initial state radiation simply adds a small, con- 
stant "background term" to the radiation pattern, which is of course independent of 
the decay width of the decaying W's. We note also that near the WW threshold there 
are kinematic constraints which are different for initial and final state radiation. The 
former is limited by the maximal kinematically allowed energy 

u < u™ ax = M w v 2 w « M w , (3.12) 

whereas the constraint on the energy emitted in the course of decay is much less severe, 

W < "ma* = ~ M W , (3.13) 

where m is the minimal invariant mass of the W decay products. 

Since we are mainly interested in the region of photon energies uo ~ V where the 
width has an important effect on the radiation pattern, we can imagine choosing a 
kinematic region, 

u ~ T < M w v 2 w < M w , (3.14) 

where the soft bremsstrahlung approximation can be used for both initial and final 
state radiation, without having to worry about the kinematic restrictions on the photon 
energy. A derivation of the radiation pattern for this situation using the classical 
picture considered in Appendix A is presented in Appendix C. 

A second difference arises when we consider extending the analysis to hadronic 
W decays, i.e. to W — > qq' . This is of course the dominant W decay channel - 
about 44% of WW pairs decay to a four-jet final state, and only about 5% of the 
decays have purely leptonic (e or fi) final states. In order to achieve a measureable 
event rate, therefore, it will probably be necessary to demand at least one hadronically 
decaying W. 

For W — > qq' both decay products can now radiate, and the pattern of radiation 
is correspondingly more complicated. Thus the electromagnetic current caused by the 
leptonic decay W + (q) — * e + (p)u e (p), 

Je ~ W) " W) ' ( } 

becomes, for the hadronic decay, W + (q) — > u{p)d{p) , 

f)^ f)^ rift 

r ^ Q W) + (1 - Q) W)-W)' (3M) 

where Q = 2/3 is the electric charge of the u quark. The general final state radiation 
pattern can then be formed from these currents in the usual way, 

A/" oc (j 1 %) 2 + (j 2 %) 2 -2xH-(ji%)(j2%) (3.17a) 
= (1 - XH) • f (Ji%) 2 + (J 2 %) 2 I + XH • U - ti) e»? • (3.17b) 
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Equivalent formulae to those considered above for gluon radiation in t decay can 
then be derived. The analysis is simplest if we assume that in each W decay the quark 
and the accompanying antiquark are anti-parallel, otherwise additional angles have to 
be introduced; this is certainly justified if the Ws are produced at rest. If we make 
the further assumption that the velocities of the quark and antiquark in the decay 
are equal, and that experimentally quark and antiquark jets cannot be distinguished, 
then we obtain the radiation patterns for the two-quark and the one-quark-one-leptonic 
decays (see Appendix D): 



(2Q-l) 2 + v 2 cos 2 



(1 + v\ cos 6*i 



+ n 2 



(2Q f - 



■l) 2 + t; 2 cos 2 9 2 

1 + f 2 COS0 2 ) 2 



+ 2 xH ■ J 



V\V2 COS#i COS 02 



Tlx 



(1 + vx cos 6*i 
(2Q-l) 2 + t; 2 cos 2 Ox 



+ 2 *H ■ J 



1 + vx cos 9i) 2 
v i cos 6i 



[1 + v 2 COS #2) 
4- ^ 2 



(3.18) 



l+Vx COS #i 



(3.19) 



where Q = Q' = | and f 1, t>2 denote the velocities of the up- type quark and antiquark 



respectively for the (qq) case, and the quantities Tlx, TI2 and J are defined in fl2.1|) and 
([2.2fl. These expressions are the analogues of the tt result ( [2.5D derived above. The 
differential cross section is obtained in the same way, with the substitution Cpa s — > a. 

To study the #12 dependence of the total photon yield (again ignoring the isolation 
cuts which will be required in practice) one has to evaluate the integrals over the 
photon radiation angle of the interference terms in (|3.18|J3.19| ). These integrals are 
finite at f ! = t> 2 = 1 and so we can use the ultra-relativistic approximation in this 
case. As shown in Appendix D, the total photon yield takes the form 



duo 



a 

7i 



LO 



'12 



(3.20) 



where we denote by (a/3) decay channels of the W + W system 



(a(3) 



(££) 
(q£) 

(qq) 



W + W 

w + w 
w + w 



ev + £is, 

lv + qq , 
qq + qq . 



The branching ratios are approximately B t 
for £ — e, fj,. We have 



(3.21) 

4/81, B ql = 24/81 and B qq = 36/81. 
In 1 " 0086,12 + 1 , (3.22a) 
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In 



sin 



12 



ln^ + 1 



with the ^-independent contributions 



mdcp. 



Z( V ) 



ilni— 

V 1 — V 



- 2 » In 



rc£L = g 2 x^) + (i-g) 2 x(t;)-2g(i-g) + 0(1-^,1-^) 



(3.22b) 
(3.22c) 

(3.23a) 
. (3.23b) 



In ( |3.23a| ), v is the velocity of a charged lepton; v and u in ( |3.23b[ ) are the velocities 
of the up-type quark (antiquark) and, respectively, down-type antiquark (quark) orig- 
inating from W + {W~). Once again, the uj dependent terms in (|3.22a|) can be either 
positive (larger #12) or negative (smaller #12). For the double-leptonic channel the 
critical angle is the same as for the it case (~ 75°) whereas with at least one hadronic 
decay channel it follows from ( |3.22b| ) and p.22c| ) that the corresponding critical angle 
is approximately 47°. 

Unfortunately, the pure "partonic" prediction ( |3.23b[ ) for the independent photon 
radiation off the quark- antiquark antenna is too naive in practice, since quark ve- 
locities are not well-defined for light quarks, and the result takes no account of the 
hadronization process where integer charge hadrons are formed and indirect photons 
from hadron decays appear. 

However these complications do not affect the main physical property of fl3.22j ), 
namely the fact that the ^-dependent (and thus the T-dependent) part of the photon 
radiation pattern is under control. As long as W + and W~ initiated jets evolve 
independently from one another, the extra yield of indirect photons remains insensitive 
to the event geometry. 

Finally, we note that an exactly analogous study could be performed for the soft 
photon radiation pattern in the process e + e~ — > Z°Z° — > fff'f. The four charged 
particles in the final state give rise to a rich interference structure, as for the four-jet 
decays of the WW pair discussed above. The four charged leptonic decays of the 
Z°Z° pair would provide a particularly clean environment in which to study this, but 
in practice the event rates would be prohibitively low. 



4 Numerical results for top 
4.1 Preliminary remarks 

In this Section we illustrate the behavior of soft radiation near threshold with some 
examples. We will discuss the gluon distributions in tt events case in some detail, 



8 



and show a single example for photons in the W + W~ case. Our emphasis will be on 
the influence of the width T on the radiation pattern. As discussed above, one might 
expect that, because of the top quark's large width (and hence short lifetime), the bb 
pair would radiate coherently, as if they had been produced directly. We will see that 
the correct soft gluon distributions can differ considerably from those which arise from 
that expectation. We will explore the sensitivity of the radiation patterns to T and 
at the end we will consider briefly to what extent this sensitivity might be useful for 
measuring T. 

Before presenting our numerical results it is helpful to summarize the results from 
Section 2 that are most relevant to what follows. Recall that the gluon emission 
probability density is given by Eq. (Ol) with 



_ v 2 sin 2 6>i v 2 sin 2 6> 2 . , t; 2 (cos6>i cos6> 2 - cos6>i 2 ) 

(l-^cosflx) 2 (l-wcos^) 2 XM (l-vcos^Xl-wcos^)' [ ' 

and x( u ) = r 2 +L 2 '■> — X — 1- The profile function \ determines the amount of 
interference between radiation off the b and b: for T = 0, we have \ = and the 
interference is suppressed, whereas for r ^> u, x ~ > 1 an d we have the full coherent 
emission mentioned in the previous paragraph. Finally, note again that 

• The only dependence on the relative geometry of the b and b appears in the 
interference term; independent emission has no #i 2 dependence. 

• In contrast to the high energy case [2] where we saw mostly destructive inter- 
ference, the interference term here can have either sign. We will see explicitly 
below that, as pointed out earlier, it tends to be constructive for small # 12 and 
destructive for large 6*i 2 . 

• From the form of x h is clear that we have maximum sensitivity to the width 
when r ~ uj (x not near or 1), i.e., when the energy of the radiation is 
comparable to the width. 

The remainder of Section 4 amounts to an elaboration of these points. 

Now, these interference width effects influence the radiation pattern, as illustrated 
in the differential distributions we present below. The behavior discussed above would 
be clearly evident if the differential distributions were observable, i.e. if we had access 
to arbitrary gluon energies and arbitrarily large numbers of events, and we could ob- 
serve partons directly. Of course, in the real world there are jets and limited statistics. 
Therefore in what follows we make some concessions to reality by also considering 
integrated distributions. Furthermore we assume that the gluon will be 'detected' as 
a soft jet, and so we take 5 GeV as the minimum observable gluon energy. However, 
this is no substitute for realistic simulations, nor is it meant to be; our results are 
meant to suggest the kinds of distributions that would be interesting to study. 
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The results presented below are more or less independent of the top mass. If 
we treat T as a parameter, m t comes into ( 4.24| ) only weakly through the b quark 
velocity v. For definiteness, we use v = 0.9944, which corresponds to m t = 140 GeV 
for m& = 5 GeV and my/ = 80 GeV. A further generalization is possible if to is 
kept fixed: then V and u only enter through \. Now, what values of x are relevant for 
accessible gluon energies and interesting top masses? For our canonical mt = 140 GeV, 
T w 0.7 GeV in the Standard Model. Taking u = 5 GeV, we obtain \ = 0.02, 
which is rather close to zero. The gluon distribution is therefore close to that for 
independent emission, not coherent emission as we might naively guess. The story is 
slightly different for larger m t because the Standard Model width grows as mf . For 
example, for m t = 200 GeV, T approaches 3 GeV and with u = 5 GeV we have 
X ~ 0.3. So the heavier the top quark, the larger are the relevant values of T and \- 



4.2 Soft gluon distributions 
4.2.1 Radiation out of the bb plane 

M takes a particularly simple form for gluons radiated perpendicular to the plane of 
the bb pair. In this case, Q\ = $2 = vr/2 and 

oc 1 — x cos Q\2- (4.25) 

The width and energy dependence are given, via Xi D Y the extent of the deviation 
from constant behavior of the gluon emission probability as a function of 9\2- This 
is illustrated in Fig. 1, where we showQ uj(dN / duodVt) as a function of # 12 for x = 
0, |, |, and 1. Interference gives rise to a dramatic difference between the independent 
emission (x = 0) and coherent emission (x = 1) cases. Note that, as indicated above, 
the distribution for x — 0.02 for our canonical 140 GeV top will be very different from 
the expected coherent emission case, and a heavier top with a larger width would lie 
somewhere in the middle. 

We can turn the discussion around to ask under what circumstances we are sensitive 
to the exact value of T. As we have emphasized above and as is clear from Fig. 1, 
we have maximum sensitivity for x n °t too close to or 1. So for m t = 140 GeV 
distributions of 5 GeV gluons are not very sensitive to a Standard Model T, but if we 
could observe 1 GeV gluons or if the width were much larger than the Standard Model 
value, x would be in a more sensitive range. For a heavy top, though, soft gluons are 
sensitive to the Standard Model width, which follows from the value of x i n the 200 
GeV example above. 

Now this out-of-plane distribution has clear, simple properties, but it is a differ- 
ential distribution - a snapshot. Let us consider to what extent the characteristics of 

-'^The factor of w in the normalization is chosen so that the only u> dependence is in x, cf. Eq. 
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Fig. 1 are retained in integrated distributions. In Figure 2 we show the distribution of 
gluon radiation out of the bb plane, now integrated over ±7r/8 in each angular direction 
and over gluon energies from 5 to 10 GeV. Because we have integrated over uj we now 
use T rather than x to label the curves. We see in Fig. 2 the same structure as in 
Fig. 1, and we may draw similar conclusions. We also note that integration over these 
angles represents ~ 20% of the gluon events. 

4.2.2 "Angular ordering" effects 

As is well-known, gluon emission off colourless qq pairs exhibits so-called angular 
ordering behavior (see for example Refs. [14,15] and the discussion above in Section 
2): If we split the radiation into pieces associated with the quark or the antiquark, 
then integrate over the azimuthal angle about, say, the quark's direction of motion, 
the quark piece of the radiation vanishes for polar angles greater than the q-q angle. 
(In particular this implies that all radiation is suppressed for collinear q and q.) This 
suppression is due to interference between gluons radiated by the q and the q. A 
derivation of this result is given in Appendix E. 

Because the width controls the interference in the case of top, we expect to see 
angular ordering behavior (or not) according to the size of T. To explore angular 
ordering for radiation off fo's from top decay, we will examine the gluon emission 
probability integrated over the azimuthal angle with respect to the b direction for 
fixed values of 9± 2 . The results for 12 = 5°, 30°, 90° and 150° are shown in Figures 
3(a), 3(b), 3(c), and 3(d), respectively. 

Before considering width effects, let us discuss the general features of these curves; 
for this purpose Fig. 3(c) is the most instructive. The distribution is shown as a 
function of 9, the gluon's polar angle with respect to the b direction. First we see, near 
9 = 0, the dead cone characteristic of emission off heavy quarks: emission is suppressed 
along the b direction but peaks nearby at the 'dead-cone' angle, 9 C ~ m b /E b m 6° for 
m t = 140 GeV. There is something like a dead cone at 9 — 9 12 = 90°, due to the b. 
The reason for the asymmetry between the b and b is that we have integrated about 
the b azimuthal angle, but we did not separate out the b radiation. A final general 
feature we note is that most of the radiation is in the vicinity of the quarks, so the 
distribution falls off at very large 9\ 2 . 

Now we examine how the width affects angular ordering. In Fig. 3(a), 9 12 = 5°. 
The azimuthally integrated distribution is shown for x — (solid line), 0.5 (dotted 
line), and 1 (dashed line). With the b and b so close together, there are not two 
distinguishable dead cones but one broadened peak. As x increases from we see 
suppression of the radiation by the interference, and for x — 1 the emission is nearly 
eliminated; recall that for collinear q and q we expect no radiation at all. As large 
as this effect is, however, it is of academic interest only: because b jets have a finite 
angular spread we can never hope to identify events with this configuration. 
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In Fig. 3(b) we consider a larger bb angular separation, 9 12 = 30°; we can just 
begin to discern the effect of the b dead cone. We see again that as x increases from 
and the interference turns on, emission at angles larger than 12 is suppressed. We get 
maximum suppression — as much as an order of magnitude — for the coherent case, 
X — 1. Between the b and 5, that is, for 9 < 9 12 , the width makes no visible difference, 
but at larger angles the interference is destructive. 

In Figure 3(c) we show the same distribution for 9\2 = 90°. The effect appears less 
dramatic: for any given 9\ the difference between x = an d X = 1 is n °t very large. 
However, on closer inspection we notice that the curves cross at 9 — 9\i- The radiation 
outside the b and b (i.e. 9 > #i 2 ) is again suppressed, but now the radiation between 
the b and b is enhanced, so that the suppression of radiation outside the bb pair relative 
to the radiation between is larger than it appears at first. The net interference for this 
larger 9 U is constructive. 

We consider a nearly back-to-back bb pair in Fig. 3(d) where 9 12 = 150°. For 
such a large angular separation there is little room outside the bb pair and the entire 
angular ordering effect amounts to an enhancement of radiation between the b and b. 
This is the well-known 'string' [18] or 'drag' [19] effect. 

As an aside, we recall that all of the effects we discuss are also relevant to photon 
radiation in the W + W~ case. A more complete numerical treatment of this will 
be given elsewhere [13], but in the meantime we give an angular ordering example 
for illustration. We show in Figure 4 the W + W~ analogues of Figs. 3(b) and (d), 
for x — (solid lines) and 1 (dashed lines). It should be noted that we have not 
included the effects of initial state radiation, so that this figure corresponds, e.g., to 
77 — > W + W~ rather than to e + e~ — > W + W~. Now because the lepton is nearly 
massless, the radiation peaks are much sharper than for the fo's, but otherwise we see 
similar features: suppression outside and enhancement between the l~ and l + for 9 U 
large and small, respectively. However, there is one other important difference: because 
of the large W width (2 GeV) and lower accessible photon energies, the relevant values 
of x reach the width-sensitive range. Furthermore, leptons do not hadronize, so small 
lepton-antilepton angles are accessible. Given sufficient event numbers, these effects 
should be observable. 

Returning to top, we have seen in Fig. 3, then, that the overall effect of the width 
on these azimuthally integrated distributions is enhancement between the b and b and 
suppression outside them. This suggests another way to get at the width dependence 
experimentally: look at the net radiation between and/or outside the bb pair. We 
show in Figures 5(a) and 5(b) distributions integrated over as in Fig. 3, over gluon 
energies (again from 5 to 10 GeV) and also over 9. In Fig. 5(a) we integrate 9 from 
to 9 12 - for the radiation between the b and b - and in Fig. 5(b) from #12 to n - for 
the radiation outside. 

In Fig. 5(a) we see the enhancement of radiation between the b and b as the width 
increases, and, as suggested in Fig. 3, the effect increases with increasing 12 . The 
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radiation outside the bb pair shown in Fig. 5(b) exhibits suppression of emission as 
the width increases. The largest effect here is at small 8u; however we cannot hope to 
do a measurement at very small angular separations because, again, b jets have finite 
angular size. In both of these figures we see, as we did in Fig. 2, that the Standard 
Model width for m t = 140 GeV, 0.7 GeV (dotted line) gives distributions very close 
to those for independent emission, and that sensitivity to T is only obtained in the 
few GeV range. 

4.2.3 Integrated distribution 

Finally we integrate over all angles and gluon energies from 5 to 10 GeV to show the 
total soft gluon emission probability as a function of #12 in Figure 6. As we shall see 
below when we discuss event rates, this may be the only distribution we have very 
much hope of seeing without multi-decades of collider runs. The independent emis- 
sion case, T = 0, is, as always, completely independent of 6> 12 due to the absence of 
interference. Increasing the width turns on the interference and induces #12 depen- 
dence. The interference is destructive at small 9 U and constructive at large 6> 12 , and 
we see the crossover point crit 75° discussed in Appendix B. Also evident here is 
the complete suppression of all radiation for the coherent case (large T) when the b 
and b are collinear. And, again, there is sensitivity to T as it approaches the few GeV 
range, i.e. as it becomes comparable to the gluon energies we consider. Finally, note 
that the gluon emission probabilities are of the order of 20%. 

4.3 Event rates 

This brings us to a discussion of event rates and the prospects for measuring the top 
width from soft gluon distributions. We must re-emphasize that our results are at 
the parton level, and any realistic assessment must incorporate hadronization of the 
6's and gluons as well as detector resolutions and acceptances, etc. Having said that, 
we now look at cross sections. The cross section for ti production near threshold for 
m t = 140 GeV is about 1 pb [21]. If we assume a yearly luminosity of 10 fb _1 , this 
implies 10 4 ti events per year. We saw in the previous subsection that the number of 
events with a soft gluon is roughly 20% of the lowest order rate, or 2000 events. If 
we further require leptonically (e or /i) decaying Ws, we are reduced to about 100 
events/year. This is less than promising, but not quite hopeless. If there is no other 
viable way to measure the top width, soft gluon distributions may be an option. 

5 Conclusions 

How well does measuring the top width from soft gluon radiation in e + e~ — > ti com- 
pare with the standard technique of scanning the threshold structure of the total cross 
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section? Each method has its disadvantages. The threshold structure is subject to un- 
certainties from beamstrahlung and beam energy spread, and from theoretical higher 
order corrections and dependence on parameters like m t and a s . The soft gluon ra- 
diation method avoids these problems, but it is a higher order process with a lower 
event rate. The two methods could therefore be considered as complementary - the 
threshold cross section loses sensitivity with increasing width, but as we have seen, 
the gluon radiation pattern becomes more sensitive at larger T for accessible gluon 
energies. For most of the expected top mass range, the threshold structure method is 
probably better, but if m t and T are large, then examining soft gluons may be more 
useful. 

In summary, we have seen that the top quark's large width gives rise to new effects 
from the interplay between the strong and weak interactions, and that the top width 
affects the distributions of soft gluons radiated in top events. Near the ti threshold, 
the effect of the width is to suppress the interference between gluons radiated by the 
b and b, in contrast to the expectation of coherent radiation from the bb pair. The 
sensitivity of the gluon distribution to T is largest for gluons with energy uj ~ T. If 
5 GeV corresponds to a realistic minimum energy for measurable gluon jets, then the 
Standard Model width of a 140 GeV top quark is sufficiently smaller than this to 
almost completely suppress the interference. 

Note that the results of our analysis could in principle be incorporated into a Monte 
Carlo scheme for generating final states in e + e~ — > tt — > W + bW~b events. Contrary to 
the standard expectation (see for example [22]), the bb antenna is practically inactive 
here since the bulk of the radiation, that is primary gluons with uj > T (in the top 
rest frame), is governed by the tb and ib antennae, and is thus unaffected by the 
relative bb orientation angle 9 bb . When parton cascades are included in the picture, 
the corresponding hard scale Q is given by 

Q ~ E b » M2 ~m™ ■ ( 5 - 26 ) 

The only particles which are sensitive to 9 bb are those originating from primary brems- 
strahlung gluons with uj < T, whose yield is determined by the parton cascade scale 

Q ~r. 

Finally, we have also extended our analysis to the case of photon radiation near 
threshold in e + e~ — > W + W~, including both hadronic and leptonic W decays. This is 
especially relevant for LEP200, where the measurement of soft photons with uj ~ Tyy 
would also reveal interesting interference effects. We will present numerical results for 
this elsewhere [13]. 
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A Semi-classical derivation of the radiation pat- 
tern 



The decay of a heavy t quark at rest produces a fast-moving b quark and thus causes 
acceleration of the colour charge. We are interested in gluon bremsstrahlung induced 
by this acceleration. Analogously to the treatment of classical electromagnetic currents 
[20] , the colour field formation is conveniently described in terms of Lienard- Wiechert 
potentials. Thus the two quark currents which participate in the colour field formation 



are: 



£ = tfi<f 3 (f-iTi(t-toi))-0(*-*oi), 
32 = v 2 5\r-v 2 {t-t Q2 ))-d{t-t Q2 ) , 



(A.la) 
(A.lb) 



where Uo are the times of the two decays. The emission amplitude for the field com- 
ponent with the 4- momentum (u, k) is proportional to the Fourier transformed total 
current, which we write introducing the effective "colour charge" as 



]{k) = ^glC~ F (Uk)-Uk) 



where the relative minus sign reflects the opposite charges, and g 
of the two terms of ( |A.2| ) we have 



(A.2) 

Ana a. For each 



i 



3i(t,r) 

ik°t 0i 



k° - (k ■ Vi) 

The field potential induced by the current (|A.3|) at large (positive) time x° reads 



(A.3) 



AAx) 



d A k 



-2m8{k 2 )]-l{k) 



d 3 k 



(2tt) 4 l J 2o;(27r) 3 

The momentum Fourier component of the total vector field is given by 

A(k) = JgJc~ F (A^k) - A 2 (k) 



(A.5a) 



Ai{k) 



Vi 



u(l — Vj cos OA 



o iu>t 0i 



uj = \k\ , {k ■ Vi) = uJVi cosOi . (A. 5b) 



To calculate the radiation probability we square the projections of the full field am- 
plitude ( |A.5a| ) onto two "physical" gluon states e\, where (e\ ■ k) =0, (ex) 2 — 1, and 
sum over polarizations to obtain 



dN 



d 3 k 



2cj(2tt) 3 



A=l,2 



d 3 k 



2w(2vr) 3 



Sat 



k a kp 
P 



A*?(k) . 
(A.6) 
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Substituting for the vector field from ([A. 5 ) and making use of the relations 



Via 



k a k/3 

Oafi — 
8a/3 



k 2 
k a kp 

k 2 



= v i sin 9i 



V2/3 = V1V2 (cos 6*12 — COS 6*i COS 6*2 ) 



(A.7a) 
(A.7b) 



we conclude that the |^4i| 2 and |^2| 2 terms reproduce the sum of the "independent" 
radiation contributions 7Z indcp . fl2.1|) while the interference 2Ke(Ai A£) is proportional 
to 2 J (|2.2|) . So finally we arrive at 



dN = — ^ • \K X + K 2 + 2 Re e M*a-*») . 



(A.8) 



This expression describes the radiation accompanying the process with heavy top 
quarks decaying at times t 0i after the tt production. These times are not measured, 
but are distributed according to the decay exponentials 



■00 

T / dto e~ 
'0 



rt 



(A.9) 



Substituting (|A.8|) into the decay-time integrals we see that the interference term gives 
the x factor, 



r / dt 0i e~ 
Jo 

r r 



r =1= iuj 



1/2 



r + %oj r — iuj 



(A. 10a) 
(A.lOb) 



leading to the final expression which is identical to the representation 

Thus we conclude that the cu-dependence of the soft radiation is due to incoherence 
induced by the uncertainty At ~ in the acceleration times of the two (6-quark) 
charges. Such a delay can be resolved by a gluon with a small wavelength 



A ~ < (At ) ~ r 



-1 



in which case (u > T) the coherence gets lost and the radiation pattern reduces to the 

sum of the two independent b and b contributions, i.e. tb and bi antennas, (see ( |2.7a| )). 
On the other hand, for wavelengths large compared to r _1 , the ( |2.7b| ) regime, the time 
delay does not affect the radiation: coherence remains undisturbed and the pattern is 
given by the bb antenna describing the point-like production of the "light" quark pair, 
just as if there was no ti stage at all. The two situations can be represented pictorially 
as shown below. 
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Oti 



Oti 



Note that the result ( |A.10b| ) for the decay profile function can easily be generalized 
to the case of different decay widths. This arises, for example, in the case of the 
production and decay of a pair of different supersymmetric particles, for example qg, 
where the decay width of the decaying particles could in principle be very different. 
In the general case we obtain 

o f r i T * \ 1^(1^2 + a; 2 ) 

xM = Re IfTT^ f^i = (H + ^)(r| + ^) • (A ' n) 

For example, for very different decay times, say Fx ^> T 2 (and ou ~ T 2 ) the expression 
( [A.llj) would lead us back to the original 

However we note, in this context, that our perturbative treatment of gluon emission 
only makes sense if u > fi ~ 1 fm~ . Unfortunately this means we cannot discuss in 
this way the interesting cases of charged Higgs decay H + — > tb or single top production 
gW + — > tb where Tb -C /i and x ^ 1- 



18 



B Angular-integrated distributions in e + e —> tt 

Using the result 



1 , v 2 sin 2 9 

(1 — V COS9) z 



In 4 



-i ( i — reoswj- v 1 — v 

we see that the angular-integrated contributions of the independent terms are 



T? 

. 'Mndcp. 
47T 



1 1 1+Vi 9 - Tf 

— In 2 = l{Vi 

Vi 1 - Vt 



The coherent contribution (|2.2|) contains the integral 



dQ 



2(1 — V\V 2 cos 9 



12) 



4tt (1 — v\ cos#i)(l — V2 cos ^2) •/ 27r (pik) (p2k) r 1 — r 



dQ uj 2 (pip 2 ) 



1,2). 



1 . 1 + r 
in 



(B.l) 



(B.2) 



where 



r = 



;i-^)(i-^l) 



M 2 Mf 



\ (1 - ^i^cos^s) 2 \ (pip 2 y 



(B.3) 



(B.4) 



the Lorentz-invariant quantity that is closely related to the relative quark velocity in 
the rest frame of the pair. In terms of the invariant energy s — (pi + p 2 ) 2 and the 
c.m.s. momentum 

[s - (M! + M 2 ) 2 ] [s - (M!-M 2 ) 2 ] 



one has 



4s 

2p c y/s 



;b.s) 



B.6) 



s - (M 2 + Mf) ' 

For the case of equal masses, Mi = M 2 , ( |B.4j ) is related to the quark c.m.s. velocity v c 

(B.7) 



by 



1 + v 2 



Combining ( |B.3| ) with the two remaining terms of (|2.2j ) which give a constant subtrac- 
tion, we obtain (cf. ( |B.2|) ) 



^ 1 , 1+ r ^ * 

T" Koher. = - In 2 = T(r 

47r r 1 — r 



;b.s) 



Putting everything together gives, for the angular integrated radiation yield, 

C F a s 



dN 
du 



TTUJ 



(i-xH) 



— log 1 log 4 

Vi 1 — Vi v 2 l — v 2 



1 1 + r 

-log- 2 

r 1 — r 



(B.9) 
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We see that the total radiation splits into incoherent and coherent contributions, 



^^{(1-XM).[IW+1W] + X (u)-l(r)}, (B.10) 

the relative weight of which is controlled by the profile function x which depends on the 
uo/T ratio. The first contribution consists of two 1(vi) terms describing independent 
radiation off quark antennae "attached" to the c.m.s. (i.e. in the rest frame of the 
decaying top quarks). The function I can be expressed in terms of the "4-angle" rj of 
the quark momentum as 

^'•^(dkr 1 ) ; « =*■*<*> ■ 

The coherent contribution, in contrast, carries no information about the initial it 
system but depends exclusively on the relative motion of the two final colour charges, 
the (12) antenna. The argument of the corresponding X factor here is the relative 
"4-angle" between the quarks, 

I M = 2 (^A)" 1 ) ; r = ta ' lh < A >- < R12 » 

where A = rji + 772 has to be calculated in a reference frame where qq are anti-collinear 
(e.g. in the c.m.s. of the pair). 

In practice, we are usually working in the ultra-relativistic limit (1 — Vi) <C 1, 
where Eq. ( P-9|) is dominated by the logarithmic collinear singularities. Since the 
angular integral of the interference term J (Eq. ( |2.2| )) converges at V\ — v 2 — 1, one 
would expect the main collinear contributions oc log(l — V{) to be tu-independent. To 
verify this let us keep only the non- vanishing logarithmic and constant terms in (|B.9| ), 
neglecting powers of (1 — Vi) <C 1. We can approximate ( P-4|) as 

9 4 (1- vx)(l -v 2 ) 

r 2 « 1 i ^ ^ ? . B.13 

[ 1 - cos 9 12 + [(1 - Vl ) + (l-v 2 )] cos 9 12 f 

Thus for #12 values not particularly close to 0, so that we can neglect the second term 
in the denominator compared to the first, we obtain 

1-^2. ' 1 -'"' (1 "? ) « 1 - < R14 > 
(1 - cos6>i 2 ) 2 



and the expression in curly brackets in (p.9|) becomes 



{} 



2 2 

In h In 4 

l-Vi l-v 2 



2x(w) 



1 - COS 012 

In h 1 



(B.15) 
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Thus we have verified that the main "collinear" contributions are ^-independent in 
the region of relative quark angles #12 exceeding the aperture of the corresponding 
"dead cones" . It is interesting to compare the relative size of the two terms in (|B.15|) 
in practice. For example, for the case of a 140 GeV top quark pair we find 



{ } ~ 7.9 + 2 X H 



1 - COS 6>i2 , . 

m h 1 



showing that the integrated quantity is indeed quite sensitive to \ and hence to the 
width. (In contrast, for a WW pair decaying to muons and neutrinos the numerical 
term from the logarithms is 23, and the sensitivity to T is much decreased.) One could 
imagine, for example, measuring the profile function x by studying the 0i2 variation 
of the total radiation yield. 

Note in particular that the second cu-dependent term in ( B.15 ) enhances the radi- 



ation at large 0i2 and acts destructively when 12 is chosen below the value G crit given 
by 

cos6 crit = l-2exp(-l) = 0.26424, 

6 crit « 75°. (B.16) 

It is interesting to notice that the suppression at 0i 2 — ► can be strong enough 
to completely compensate the main collinear contributions. Indeed, taking paramet- 
rically small angles 

0? 2 « (1- Vl ) + (l-v 2 ), 
we would obtain for r the value (see ( |B.13| )) 

r 2 « 5 (B.17) 

[(l-^i) + (l-^ 2 )] 2 

which can be arbitrarily small for nearly equal quark velocities. If r <C 1 the sec- 
ond (coherent) term in the general expression (|B.9|) for the radiation yield becomes 
negligible and one is left with (v = v\ ~f 2 ) 



dN C F a s 

— » 2 1-xH 

d(jJ 7TU1 



ih-H— 2 

V 1 — v 



The result vanishes when xi 00 ) ~ * 1 (*- e - f° r 00 ^ T). This corresponds to the total 
coherent suppression of radiation off two opposite charges moving in the same direction 
with equal velocities. Notice that for uj = T the independent radiation off collinear 
daughter particles is suppressed by a factor 2. 

It is worth mentioning that the tu-dependent coherent effect cancels after inte- 
gration over all angles 8u. This means that the interference J does not affect the 
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total bremsstrahlung caused by the decay of the heavy unstable objects, but only re- 
distributes the accompanying radiation between configurations with different relative 
angles $12 . This can be checked explicitly by evaluating the angular integral of the 
second (coherent) term of ( |B.9j ) over 12 which results in 



sin 9i 2 d6i2 
2 ■ 



dQ 

47T 



dcos 612 



ihl±I-2 

r 1 — r 



— In — — — - H In 

Vi 1 - v 1 v 2 



1 , 1 + v 2 A r . n in f dfl _ 

- m 4 r = / S111 ^12^12 ■ / ^indep 

) 2 1 — V 2 ) JO J 47T 



(B.18) 



This suggests a way of extracting the profile function by studying the coherent 
"redistribution" effects in the total radiation yield. This could be done for example 
by comparing the gluon yield at fixed #12 below and above the crit value. 

Let us consider, therefore, the integrated quantity characterising the interference 
effects, namely the difference of the integrals over the bb opening angles above and 
below crit : 



r _ 1 r ■ n M ^ dN 1 

6 = TT< / sm9 12 d9 12 — — 

(1 + cos H crit ) Je crit du (1 — cos H c 



e cr it ujdN 

sin 12 d0 12 (B.19) 
duo 



The normalization in (p.!9|) is chosen so that the 812-independent contributions cancel. 
Within the relativistic approximation, we have for ( |B.15| ) 



M = n indcp . + 2 X (co) 



1 — COS 019 

In + 1 



(B.20) 



Making use of the definition of the critical angle ( |B.16| ), we calculate the integrals 



cos e crit 



^1 + COS crit ) J-l 
1 rl 



(1 — COS© crit ) JcOsScrit 

to obtain finally for ( B.19|) 



d cos 9i2 
d cos 9i2 

C F a s 



1 - COS 6»i2 

In h 1 



1 - COS 012 

In h 1 



1 - cos Q crit 
1 + cos crit 



7T 1 + COS 6 crit 

Invoking the numerical value of the critical angle ( |B.16| ) we write ( |B.22| ) as 

a s (uj) 



4.218 



7T 



,(B.21a) 



(B.21b) 



(B.22) 



(B.23) 



Notice that we have chosen here the gluon energy as the argument of the running 
coupling, since it is large gluon energies that contribute to ( |B.23| ). 
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C Initial and final state radiation in e + e 
near threshold 

Let us denote by vq ( — ^b) the 3-velocities of the incoming e + (e - ). 
state electromagnetic current takes the form 



w+w~ 



Then the initial 



Jo 



v S 3 (f-v t) - {-v )S 3 {f +v t) -&(-t) 



which gives an extra contribution to the field amplitude (see Appendix A) 



An 



+ 



1 - (v n) 1 + (v n) 



v 



UJ 1 — Vq COS 2 6q 



(C.l) 



(C.2) 



where n is direction of the photon and 8q its angle with respect to the incoming 
positron. Since the e + e~ "disappear" at the same time as the WW pair is produced, 
the amplitude ( |C.2| ) is real with our convention (t — 0). Let us choose for the sake 
of simplicity the {Ivlv) decay channel of the WW. We use ( |A.6| ) to obtain for the 
radiation probability 



dN 



d 3 k 



2w(27r) 3 
a(k) = — 



E \Mk)-e x 



A=l,2 



■ e 



v 2 



Win) 



1 



• e 



0,(9=1,3 



2v 



[v 2 n 



[v n) 



(C.3) 



where V\, v% stand for positively and negatively charged final state leptons respectively. 
The novel feature of (p.3| ) is an interference between the initial state radiation (ISR) 
(the last term) and the final state radiation (FSR) (the first two terms), which has 
the structure 

Joi ■ 2 Re e iuJt(n +J 02 -2 Re e iuJt ° 2 . (C.4) 

After the integration over the W decay times £oi is performed (see ( |A. 1C ) ) it gives 
rise to the same profile function x, due to the identity 



Re 



r 



Re 



r 2 



r + iuj r — %u J l r ± %u J r 2 + lJ 1 

The final answer can therefore be represented as 

dN = M-s+Mb+M/f 

where the first term is, as before (cf. (|2.5|)), 

A/" FS 



n 1 + 'Ji2 + 2 X {uo) 



(vin)(v 2 n) - {vA) 
[l-(^n)][l-(^ 2 n)] 



(C.5) 



(C.6) 



(C.7a) 
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The new terms describe independent radiation off the initial state e + e antenna, 

jv- = n v • < a7b > 

(1 — Vq cos 2 9 ) 2 
and the ISR/FSR interference contribution 

1 - {v n) z { 1 - (vin) 1 - (v 2 n) J 

Note that the ISR/FSR interference (|C7c[ ) vanishes (i) after integration over the angles 
between the ISR and FSR antennae (keeping the relative angle between the daughter 
charged particles fixed), and (ii) in the limit u ^> T, as expected. We recall also that 
both J\f IS and M I/F vanish when the kinematic limit uj > u)™^. (Eq. ( |3.12|) ) is exceeded. 

In practice we are not interested in photon emission close to the beam direction, 
and so for # 3> m e /Mw we can set vq = 1 and the above expressions become 



and 



Ms * -^r> (c.8) 

sin On 



\r ( \ \r f CQS ^o-(^i) cosgp - (v v 2 ) \ 

A/j/F = ' Mjb ■ < 77^7 7T^— > . C.9 

1 - (vin) 1 - (v 2 n) 
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D Hadronic W + W decays 



In this Appendix we derive the expressions for the photon radiation pattern in e + e~ — > 
W + W~ when at least one of the W's decays hadronically , i.e. to qq'. Consider first 
the case when both W's decay hadronically. The two currents are then (see Eq.( |3.16| )) 



J i 



ft 



Q 



Q' 



A 1 
Pi 

(kpi. 
n 

P2 

(kp 2 



+ {i-Q) 
+ (i - Of) 



Pi 



(kpi) (kg) 



-t 1 



(kp 2 ) (kq) 



(D.l) 



where Q = Q' = | and pi, P2 denote the momenta of the up-type quark and antiquark 
respectively. In the 3- vector form ((eg) = 0) we can write 

Q l-Q 



J 



v 



nv 



nv 



v (2Q - 1) + nv 
uj 1 — (nv) 2 



(D.2) 



where we have treated the quark and the accompanying antiquark momenta as anti- 
parallel and of equal mass. Then 



Af {qq) 



oc 



[evx 



(2Q - 1) + nv x 



1 — [nviY 
-2x(^) • (eui)(eu 2 ) 



+ 



[ev 2 ) 



(2Q' -l)+nv 2 



,nv 2 



(2Q-l)+nv 1 



[nv 1/ 



(2Q'-l)+nv 2 



[nv 2 



(D.3) 



For the mixed one-quark-one-leptonic decay configuration, Q' — 1 and we have 

(2Q-l)+nv 1 



[evt 



1 — (nvi) 2 
2x(uj) ■ (evt)(ev 2 ) 



- 2 




+ 


(ev 2 ) ^ 
1 — nv 2 . 



(2Q-l)+nv l 



[nvi 



nv 2 



(DA) 



Experimentally it seems very difficult (if at all possible) to discriminate between the 
two jets that originate from, say, W + — > u + d decay. Without being able to separate 
quark and antiquark jets (Q from 1 — Q) we have to drop in ( p.3|JD.4| ) the odd terms 
in 2Q — 1 (2Q' — 1). After summing over photon polarizations, Eq. ( |A.7|) , we arrive at 



M {qq) ocv 2 sin 2 0! 



(2Q-1) 2 + y 1 2 cos 2 1 , 2 . 2 , (2Q'-l) 2 + v 2 cos 2 9 2 



+ Vn sin 9 2 



1 — v 2 cos 2 9i) 2 

+ 2x(uj) ■ V\V 2 (cos0i COS0 2 — COS 9 12 



'1 — V 2 COS 2 9o) 2 



Vj cos6>! 
1 — v\ cos 2 9\ 



V 2 COS 2 
— V 2 COS 2 2 



(D.5a) 
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(1 - COS 2 6»i) 2 (1 - V 2 COS0 2 ) 2 /t~\ r"U N 

+ 2x(uj) ■ v\v 2 (cos 0i cos 9 2 — cos 6 ) i 2 ) 



1 — v 2 cos 2 9\ 



1 — V 2 COS 2 



Recalling the expressions for IZi, 1Z 2 and J introduced in Section 2, it is straight- 
forward to cast this result in the form given in Eqs. ( |3.18|jrT9| ). Note that there is 
a subtlety in the case concerning the definition of the angles. Whereas for the 

ti case we could unambiguously define, say, X with respect to the b quark, with two 
indistinguishable jets from the W we lose this capability. However having made a 
choice for defining X , the definition of 9 2 is correlated with that of 9\ 2 . The invariance 
of the AA") distribution under this symmetry is manifest by the dependence on the 
quadratic terms cos 2 9i, cos 2 9 2 and cos0i cos 9 2 only. 

As long as the photon direction n is kept away from the "dead cones" of the final 
charges, a simplified version of Eqs. (p.5|) can be used in which the velocities are set 
to 1: 

*/-(«) _ | + cos2 0i . | + cos 2 2 (cos 8i cos 9 2 - cos 9 12 ) cos 9 1 cos 9 2 , , 

sm 9i sin 9 2 sm 9\ sm 9 2 

\r(qt) l + cos 2 ^ sin 2 6» 2 (cos0i cos0 2 - cos0i 2 ) cos0i . 

W q ) = . 2q — + T , wri + 2xM ■ 2 a n T\ -(D.6b) 

sm 9\ (1 — cos0 2 ) 2 sm 9\ (1 — cos0 2 ) 

Expressing the relative angle 9\ 2 between the final jets (the jet and the lepton) in terms 
of photon angles with respect to the 1,2 directions, we can write the final expression 
for the angular pattern of photon radiation as 

jyfta) oc ^ + C 2 °f 9l + ^ + C 2 °f ° 2 - 2xM cot 8 X cot 9, ■ cos 12 , (D.7a) 
sin 9\ sin #2 

jyfeQ I + co f 01 + cot 2 d _l _ 2 / w ) cot ^ cotan ^ . COS 12 (D.7b) 
sin 0i 2 2 

where 0i 2 is the azimuthal angle between v\ and v 2 projections onto the plane orthog- 
onal to the photon direction, n. 

To study the 9\ 2 dependence of the total photon yield one has to evaluate the 
integrals over the gluon radiation angle of the interference terms in ( |D.5| ). These 
integrals are finite at V\ — v 2 — 1 so we can use the relativistic approximation ( p.6| ) 
for this purpose. A straightforward calculation leads then to identical results for the 
quark-quark and quark-lepton channels, namely 

dfl (cos 0i cos 02 — cos $12) cos 0i cos 02 

47T sin 2 0i sin 2 2 . . 

d£l (cos 0i cos0 2 — cos0i 2 ) cos 0i , sin0i 2 

— o ; ; = m + 1 . 

4tt sin 2 0i (1 -cos 2 ) 2 
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Taken together with the independent contributions, the total photon yield then takes 
the form given in Eqs. (gg^ggBJ). 
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E Azimuthal angle integrations 



In this Appendix we derive the expressions for the radiation pattern in tt production 
when the gluon is integrated over its azimuthal angle with respect to the direction 
of the 6-quark {i.e. the "1" direction). 

We first note that the dependence enters when we substitute 

cos 6*2 = cos 6*i cos 6*i2 + cos <j) sin 9\ sin 6*12 (E. 1) 

into the result for 7Z cohcI given in Eq. ( |2.2| ). Consider first the azimuthal average of 
the interference term, (J). We need the basic integral (a > \b\) 



2tt 



2n a + b cos <p ^Ja? — b 2 



(E.2) 



where 

a = 1 — v 2 cos Q\ cos (9 12 , b = — v 2 sin 9\ sin # 12 . (E.3) 

This gives 



V 



A 2 + sin 2 (1 - v 2 ) 2 , (E.4) 



1 — v 2 cos 9 2 , 
where 

A = cos 61 — v 2 cos 6*12 , (E-5) 

which then leads to 

(J) = ^— [y-A-cos^] . (E.6) 

1—Vx COSt^i 

As a function of 9\ the first term in the bracket is a smooth step-function-like distribu- 
tion, falling from +1 at Q\ = 0, through at cos6*i = v 2 cos 6*12 (i.e. when 9\ coincides 
with the direction of the other 6-quark), to — 1 at 9i = ir, see also Ref. [15]. 

For the azimuthal average of the TZ 2 contribution to 7Z cohcI . we can make use of the 
following: 

V = (- 1 Q , 2 ) = (v 2 ^- + 1 j V = (1 - v 2 cos6 1 cos6 12 ) ■ V 3 , (E.7) 

\(1 -v 2 cose 2 ) 2 / \ dv 2 J 

from which follows 

(n 2 ) = 2V - 1 - (1 - vl) V . (E.8) 
Note that in the special case when 8u = we obtain 

(J) = - T . V -^i ^ , (E.9) 

(1 — V\ cos ui)(± — V 2 COS &i) 
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and 



('^-'coher. ) 



(vi ~ v 2 ) 2 sm 2 e x 



(E.10) 



(1 — 1>1 cos #i) 2 (l — V 2 COS 6*i) 2 

Evidently (7Z cohcI .) vanishes when vi = v 2 , because of complete destructive interference. 

A similar analysis can be performed for the more symmetric case of the azimuthal 
average around the bisector of the "1" and "2" directions. If B is the angle between the 
gluon and the bisector and A = \Q\ 2 then a straightforward but tedious calculation 
gives for the independent contributions 



{K i ) = 2U i -l-{l-v 2 )U' i 

and for the coherent contribution 
2(1 - v x v 2 cos2A) 



(E.11) 



("^■coher. ) 



v l + v 2 ~ 2V\V 2 cos COS A 



[ v 1 U 1 + v 2 U 2 ] - (l-vl)U[ - {l-v 2 2 )U' 2 , (E.12) 



where 



Ui 

u' 



1 



1 — Vi cos 0j 

1 



;i -v 2 )sin 2 A + (cosA -^cos6) 2 J 2 , (E.13a) 
= + u i = cos9 cos A)C/f.(E.13b) 



(1 - Vi cos6j) 2 / 

It can be further shown that for the equal velocities case (v ± = v 2 = v 

2v 2 sin 2 A 



(^■coher. ) 



1 — v cos cos A 
oc A 2 , A -> 



2U~ l - sin z 6(1 - if 



U 



while for V\ ^ v 2 and 6*12 = we again arrive at (cf. Eq. ([E.lOj) 

(vi -v 2 ) 2 sin 2 6 



(1 — Vi cos 6) 2 (1 — v 2 cos 6) ; 



(E.14) 



(E.15) 
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Figure Captions 

[1] Soft gluon distribution in e + e~ — > tt — > bbWW near tt threshold, for gluons 
perpendicular to the bb plane, for x as marked. dN/doodfl is given by Eqs. fl2.4p 
and ( f4.24| ), with 61 — 62 = 90°. 6*12 is the angle between the b and b. Here and 
in all subsequent figures, unless otherwise noted, v = 0.9944, which corresponds 
to m t = 140 GeV for m b = 5 GeV and m w = 80 GeV. 

[2] Distribution of soft gluons radiated out of the bb plane, with energy and angular 
integrations, for T as marked. The distribution shown is given by 

f rlO GeV /-57T/8 /-57T/8 

dN = 2 duo sin 9 d9 d<p (dN/doodQ), 

J J 5 GeV J3it/8 J3it/8 

where 6 and are, respectively, the polar and azimuthal angles with respect to 
the b quark direction; the bb pair defines the xz plane. The factor of 2 accounts 
for radiation on both sides of the plane. V = 0.7 GeV (dotted line) corresponds 
to a Standard Model top with m t = 140 GeV. 

[3] Angular ordering effects in the soft gluon distribution with azimuthal integration: 

2tt 



dN = 00 (dN/dood cos 6) = / d<\> u{dN / dudVt). 

Jo 

The angles between the b and b are (a) #12 = 5°, (b) #12 = 30°, (c) #12 = 90°, and 
(d) 9i2 = 120°. Solid lines: \ — 0; dotted lines: x — 0-5; dashed lines: x — 1- 

[4] Analog of Fig. 3 for the W + W~ case: angular ordering effects in the soft photon 
distribution in 77 -> W + W~ -> for 6 12 = 30° and 150°. J dN = 

ui{dN/duid cos 6) as in Fig. 3. Solid lines: x — 0; dashed lines: x — 1- (The 
77 — > W + W~ case is equivalent to the e + e~ — > W + W~ case with no initial state 
radiation.) 

[5] Soft gluon distribution (a) "between" and (b) "outside" the bb pair. In (a), 

r rlO GeV i-2-k f e 12 

dN = duo d(j) sin 9d9(dN/doodQ), 

J J5 GeV JO JO 

and in (b), 

/rlO GeV r2ir pit 

dN = duo d<p sin 6d9(dN/duodQ). 

J5 GeV JO Jd 12 

Solid lines: T = 0, 00 (as marked); dotted line: T = 0.7 GeV; dashed line: 
T = 3 GeV; dot-dashed line: T = 5 GeV. 

[6] Net emission probability for gluons with energies from 5 to 10 GeV: 

//•10 GeV /■27T pit 

dN= duo d<f> sin 9d9 (dN/dudQ). 

J5 GeV JO Jo 
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